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Q . ABSTRACT 

The constructive role of non-Gaussian random fluctuations is studied in the context of the passage over the 
dichotomously switching potential barrier. Our attention focuses on the interplay of the effects of independent 
sources of fluctuations: an additive stable noise representing non-equilibrium external random force acting on 
the system and a fluctuating barrier. In particular, the influence of the structure of stable noises on the mean 
escape time and on the phenomenon of resonant activation (RA) is investigated. By use of the numerical Monte 
Carlo method it is documented that the suitable choice of the barrier switching rate and random external fields 
y ■ may produce resonant phenomenon leading to the enhancement of the kinetics and the shortest, most efficient 
reaction time. 
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1. INTRODUCTION 

3 ■ The problem of noise induced phenomena is of special interest in non-equilibrium natural systems covering broad 
class of examples ranging from physics and chemistry to biological sciences. 1 In particular, systems driven by both, 
external noises and thermal fluctuations have been shown to manifest several counterintuitive phenomena leading 
to e.g. the enhancement of the detection of weak sensory signal or optimization of a directed transport. 2 Among 
various types of "resonant behaviors" , so called resonant activation (RA) has attracted a special attention as a 
feature problem describing kinetics in fluctuating environments subject to nonequilibrium fluctuations. 3-7 The 
phenomenon is detected as a most efficient, shortest passage over a fluctuating barrier with a maximum resonant 
OA [ activation being a function of the correlation time of the barrier noise. The typical theoretical models analyzing 
the RA phenomena are usually based on a Langevin equation approach assuming the overdamped limit. 2,3,6 ' 7 
Accordingly, the influence of the external thermal bath of the surroundings on a Brownian particle is described 
in such an equation by time-dependent random force which is commonly assumed to be represented by a white 
Gaussian noise. That postulate is compatible with the assumption of a short correlation time of fluctuations, 
"j***; much shorter than the time-scale of the macroscopic motion and assumes that weak interactions with the bath 
lead to independent random variations of the parameter describing the motion. In more formal, mathematical 
Q terms Gaussianity of the state- variable fluctuations is a consequence of the Central Limit Theorem which states 
that normalized sum of independent and identically distributed random variables with finite variance converges 
to the Gaussian probability distribution. If, however, after random collisions jump lengths are ruled by broad 
distributions leading to the divergence of the second moment, the statistics of the process changes significantly. 
O ■ The existence of the limiting distribution is then guaranteed by the generalized Levy-Gnedenko 8 limit theorem. 
According to the latter, normalized sums of independent, identically distributed random variables with infinite 
variance converge in distribution to the Levy statistics. At the level of the Langevin equation, Levy noises are 
generalization of the Brownian motion and describe results of strong collisions between the test particle and the 
surrounding environment. In this sense, they lead to different models of the bath that go beyond a standard 
"close-to-equilibrium" Gaussian description. 10-12 

In the recent studies 7 we have shown that the same kind of behavior which was observed in the systems coupled 
to the sources of Gaussian fluctuations can be observed when instead of equilibrated or close to equilibrium heat 
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Figure 1. A model of the potential barrier studied in the paper. The barrier height fluctuates dichotomously between 
two barrier configurations characterized by two height H±, in this presentation H- — 0,H + — 8. A particle starts its 
diffusive motion, in t = 0, at a reflecting boundary, x — 0, and continues until the absorption at the absorbing boundary, 
x = 1, at t = t. 



reservoirs a non-equilibrated thermal baths are considered. On a single particle level, the contact with non- 
cquilibrated thermal bath can be described by a generalized Langevin equation, i.e. Langevin equation in which 
instead of Gaussian fluctuations some heavy tailed Levy-type fluctuations are incorporated. 13 Here we report 
some basic results for resonant activation processes driven by Levy flights 11 ' 12 and compare the appearance 
of the phenomenon for a Gaussian and for a heavy-tailed a-stable noise. Furthermore, we design a strategy 
that allows to determine the feature of underlying noises and recognition of their characteristic properties. In 
Section 2 the archetypal model and basic mathematical introduction are presented. More formal definitions and 
necessary technical information are moved to the Appendix. Section 3 contains results of numerical analysis. 
The discussion and final remarks are presented in Section 4. 



2. MODEL 

We study the motion of an overdamped Brownian particle moving in a potential field between reflecting x = 
and absorbing x = 1 boundaries in the presence of noise that modulates the barrier height. Time evolution of a 
state variable x(t) is given in terms of the generalized Langevin equation 

^ = —V'{x) + gt](t) + V2((t) = -V' ± {x) + V2<Z(t), (1) 

where prime means differentiation over x, ((t) is a white Levy process 8 originating from the contact with non- 
equilibrated bath and rj{t) represents a Markovian dichotomous noise of intensity g taking one of two possible 
values ±1. Both £ and rj noises are assumed to be statistically independent and autocorrelation of the dichotomous 
noise is equal to — (v))^^) ~ (v))) = ex P( — ^l\t ~ t'\). For simplicity, a particle mass, a friction coefficient 

and the Boltzman constant are all set to 1 . The time-dependent potential V± (x) is represented by a linear barrier 
with its height switching between two values H± with an average rate 7 (Fig. 1) 

V±(x)=H±x, g= H -~ H+ . (2) 



Equivalent to the stochastic ordinary differential equation (1) is a deterministic fractional Fokker-Planck 



equation 13 (FFPE) for the distribution function 



dp±(x,t) d 
dt dx 



dv±(x) , : 



DV>± {x, t) + 1Pt (x, t) - 1P± (x, t). (3) 



In the above FFPE p±(x,t) are probability density functions for finding a particle at time t in the vicinity of 
x, while potential takes the value V±(x). The fractional derivative V Q is understood in the sense of the Fourier 
Transform 14 ' 15 

/dk 
—e tkx \k\ a , (4) 

with a — 2 corresponding to the standard Brownian diffusion. The coefficient D represents the generalized 
diffusion coefficient with the dimension 13 [D] = cm a sec~ 1 . The initial and boundary conditions for Eq. (3) are 

P±(x,t)\ t =o = ^S(x), p±(x,t)\ x=1 = 0. (5) 

The first condition represents equal chances of finding initially a potential barrier in any of two possible configu- 
rations (±). The second condition is implied by the absorbing boundary, which is located at x = 1. Furthermore, 
as mentioned earlier, at x = reflecting boundary is located. 

In the approach presented herein, instead of solving Eq. (3), information on the system are drawn from the 
statistics of numerically 8, 16 generated trajectories satisfying the generalized Langevin equation (1). At a single 
trajectory level sampled from the Monte Carlo 17 study of the problem, the initial condition for Eq. (I) is 

x(t)\ t=0 =Q, (6) 

i.e. initially particle is located, with probability equal to 1, at the reflecting boundary and potential barrier is, 
with probability 1/2, in any of two possible configurations (±). Eq. (6) is a discrete analogue of the condition 
(5) for the FFPE (3). 

The quantity under consideration is the mean first passage time (MFPT), i.e. the average time that particle 
spends in the system before it becomes absorbed. The information about the MFPT is drawn from the statistics 
of numerically generated trajectories satisfying the generalized Langevin equation (1), i.e. the examined MFPT 
is estimated as a first moment of the distribution /(r) of first passage times (FPT),* obtained from the ensemble 
of simulated realizations of the stochastic process in question. More precisely, FPT is calculated as time r for 
which a Brownian particle starting initially from x = reaches x = 1 for the first time and MFPT is estimated 
as (t). 

Based solely on the statistics of generated trajectories we estimate the mean first passage time as 



MFPT=(t)=/ dx \p-(x,t)+p+(x,t)]dt= G(t)dt. (7) 
Jo Jo Jo 

where G(t) — J Q \p_ (x, t) +p + (x, t)]dx is a survival probability, i.e. the probability that at the time t a Brownian 
particle is still in the system, i.e. in the interval [0, 1). Survival probability function, G(t), is easy to construct 
from sample realizations of the process given by Eq. (1). The function G(t) could provide method for identification 
of underlying stochastic process C(t) in the cases when methods based on examination of MFPT fail to provide 
insight into underlying dynamics. Furthermore, from G(t) it is possible to calculate MFPT and to construct 
fir). 

Tn the presentation t and r are used. The r is used when MFPT = (r) is discussed or some remarks on FPTs are 
expressed. The t is used in the context of integration variable or in any context associated with the survival probability 
G(t), which is constructed by analysis of FPTs, r. This distinction is provided to distinguish time, t, which can take any 
positive values, from the FPT, r, which is a discrete random variable taking also positive values only. 



3. RESULTS 



In the standard Gaussian scenario phenomena expressing constructive role of the noise are examined as a function 
of the noise intensity a 2 , which is related to temperature. For the additive white Gaussian noise, a 2 which is 
equal to the variance of the Gaussian distribution, is the only one noise parameter. In the Gaussian regime 
stochastic resonance takes place when for a given value of a 2 examined quantity (e.g. signal to noise ratio, 
spectral amplification) is maximal. 2 For the non-Gaussian noises situation is more subtle. It is caused by the 
fact that stable distributions, which are considered here, are characterized by four parameter family (a, (3, a, /j,). 
The stability index a (a E (0, 2]), skewness parameter (3 (/? S [— 1, 1]), scaling parameter a (a £ R + ) and location 
parameter 1), The stability index describe how heavy tailed distribution is, i.e. PDFs asymptotically 

behave like The (3 parameter describe skewness of the distribution, i.e. which part of the distribution 

is "heavier", e.g. for a < 1 and (3 = 1 the distribution is totally skewed and random variable L a> i(x; a, fi) 
takes only values greater than /i. The scale parameter a scales the distribution width. All of those parameters 
determine how likely large fluctuations can be. Finally, the \i parameter is responsible for the location of the 
modal value. 

In the following simulations value of /i has been set to and a = 1/ \f2. Such a choice of a reconstruct standard 
normal distribution iV(0, 1) for a = 2, i.e. the well known Langevin equation with the additive Gaussian white 
noise is recovered and comparison with previous studies of RA 3 ' 6 - 7 can be performed. Remaining parameter, 
a, has taken values from the whole allowed range, however very small a, a < 0.2, has not been investigated. 
Furthermore, due to numerical instability a = 1 was excluded from the analysis. The MC simulation has been 
performed for a increasing by 0.1 from simulation to the simulation. Simulation for other values of parameters 
were also performed and results have been presented elsewhere. 12 

In the presented study, for a given barrier setup, various types of a-stable noises were investigated. The 
potential barrier was switching between different heights H_ = and H+ = 8. Such a choice corresponds to 
the changing of the potential barrier between the barrier and no barrier, i.e. diffusion in a potential field and 
free Brownian motion. Specific values of H± correspond to the previous study of RA 6 ' 7 ' 11 and allow for simple 
comparison with previously calculated MFPTs for the considered model in the presence of Gaussian, Cauchy 
and Levy-Smirnoff noises. 11 

For the testing purposes of the implemented numerical procedures MC solutions of Eq. (1) for a = 2 were 
compared with exact solutions of backward Fokker Planck Equation. 3, 6 The obtained numerical results are in 
the perfect agreement with the exact solutions. 3, 6 Error bars which are shown in Fig. 2 represent deviation 
from the mean and were calculated by standard MC error analysis method. 17 Numerical results presented in the 
Fig. 2 were calculated by direct integration 16 of Eq. (1), with respect to the general a-stable measure, with the 
time step dt = 10~ 4 and averaged over N — 10 3 sample realizations. The same time step, dt, and ensemble of 
N trajectories were used to analyze other cases under the study, i.e. other types of stable noises (Figs. 3-8). 
Different time steps, number of repetitions were examined, but they do not changed results quantitatively. For 
example, results calculated for dt = 10~ 5 are the same like for the larger time step, dt = 10~ 3 or dt = 10~ 4 , 
except the a = 1,(3 ^ cases, which for that reason were excluded. Presented results are divided into two 
groups. The first group correspond to the symmetric, (3 = 0, a-stable noises (Figs. 2 and 3). The second group 
to the totally skewed, (3 = 1, stable noises. 

By examination of Figs. 3 and 4 (left panel) it can be seen that the typical Gaussian behavior of MFPT 
curves changes with the change of the stability index a. For symmetric a-stable noises, (3 = (Fig. 3) with 
decreasing of the stability index, a, minima become less distinct and finally for small a phenomena of the RA 
vanishes. For totally skewed a-stable noises situation is more complicated and less regular. Obviously for a = 2 
the Gaussian limit is recovered and the RA is well visible. With decreasing stability index a minima become less 
pronounced and for a w 1.5 the RA vanishes. For a = 0.9 the RA is visible again and disappears for a < 0.7 
(right panel of Fig. 4). 

From Figs. 4-6 it can be concluded that the same kind of behavior, i.e. a typical non-monotonic shape 
of MFPT curves, is observed in two opposite limits, i.e. for equilibrium Gaussian fluctuations and for non- 
equilibrated fluctuations described by heavy tailed, asymmetric stable noise. Observation of MFPT curves only 
(Figs. 5 and 6) do not provide any information about underlying noise. Therefore more detailed study is required. 



Figure 2. MFPT(7) for linear potential barriers switching between different heights H±: '+': H- = 4, H+ — 8; 'x': 
H- — 0, H + — 8; '*': H- = —8, H + = 8. The driving noise is the Gaussian noise, which is the special case of 
the Levy stable noise, i.e. Levy stable noise with the stability index a = 2 and the scale parameter a — l/v2. Solid 
lines represent exact solution constructed by direct integration of the backward Fokker-Planck equation. 3, 6 Numerical 
results were obtained by use of Monte Carlo simulation of Eq. (1) with time step dt — 10~ 4 and averaged over iV = 10 3 
realizations. Error bars represent deviation from the mean and usually remain within the symbol size. Additionally, exact 
asymptotic values of MFPT are presented, i.e values of MFPT(7 -> 0) and MFPT(7 -> co). 
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Figure 3. MFPT(a, 7) for (3 = 0, a — l/y/2 for the linear potential barrier switching between H- = 0, H+ = 8. 
The results were calculated by direct integration of Eq. (1) with the time step dt = 1CF 4 and averaged over N — 10 3 
realizations. 
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Figure 4. MFPT(a,7) for H- = 0, #+ = 8 and /3 = l,cr = l/\/2 (left panel) and sample cross-sections MFPTX7) 
for various a: '+': a = 0.2; 'x': a = 0.8; a = 0.9; a = 1.1 and a = 1.7 (right panel). The results 

were calculated by direct integration of Eq. (1) with the time step dt — 10~ 4 and averaged over iV = 10 3 realizations. 
Additionally, asymptotic values of MFPT are presented, i.e values of MFPT(7 — » 0) and MFPT(7 — » 00). Asymptotic 
values of MFPT were calculated using Monte Carlo methods with dt = 10~ 4 and averaged over N = 5 x 10 3 realizations. 
Due to numerical problems results for a — 1 are not presented. 
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Figure 5. Sample cross section of the MFPT(a,7) surface (left panel of Fig. 4; H- = 0, H+ = 8 and (3 = l,a = l/y/%) 
with asymptotic values of MFPT, for a = 2, i.e. for the case when additive fluctuations are of the Gaussian type. 
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Figure 6. Sample cross section of the MFPT(a,7) surface (left panel of Fig. 4; H- = 0, H + = 8 and (3 = l,a = l/y/2) 
with asymptotic values of MFPT, for a = 0.9, i.e. for the case when additive fluctuations are of a heavy-tailed type. 

In Figs. 7 and 8 survival probabilities, G(j,t), are presented in a semi-log scale. Right panels present sample 
cross sections of G(j, t) surface, for various values of 7: low, high and resonant switching frequency. 

In the case when underlying fluctuations are Gaussian (Fig. 7) it can be concluded that G(7, t) for high 
values of 7 is distributed exponentially while for small values of 7 two distinct slopes are visible. These slopes 
correspond to the various time regimes responsible for passages over barrier in (±) configurations. 7 For large 
values of the 7 survival probability distribution is exponential and particle escapes through average potential 
barrier. 7 For non-Gaussian fluctuations (Fig. 8) G(7, t) is no longer distributed exponentially. Here again two 
time scales of the escape process over a fluctuating potential barrier are visible for small 7 and one for large 
values of 7. Inspection of insets demonstrates different character of G(t) for Gaussian (Fig. 7) and non-Gaussian 
statistics (Fig. 8) with noticeable difference in shape caused by different underlying fluctuations. Therefore 
methods based on the survival probability analysis can be applied for recognition of the shape of the underlying 
fluctuations. 

4. DISCUSSION 

In this study, we have investigated motion of a Brownian particle in a potential field subjected to Markovian 
dichotomous fluctuations in the presence of external, heavy tailed fluctuations arising from the contact with a 
non-equilibrated thermal bath. The potential field has been approximated by a linear potential barrier. For 
such a model phenomena of resonant activation have been investigated. Furthermore we focused on method of 
recognition of underlying noise type. 

For the numerical study we applied numerical procedures that allow to generate family of strictly stable Levy 
variables 8 and in consequence allow to perform integration of stochastic differential equations with respect to the 
family of strictly stable measures. Using that procedures resonant activation and survival probability for various 
types of stable noises were examined. In particular it was checked that for a — 2, the well studied Gaussian case 
is recovered. Monte Carlo results for Gaussian noises constructed by general procedure, for simulation of the 
generalized Langevin equation are in perfect agreement with previous numerical 6 ' 7 and analytical results. 3 For 
easier comparison with results of previous studies value of a was set to 1/ \f2. 

For nonsymmetric stable noises, for H— = 0, H+ = 8, for very small values of the stability index a RA is not 
visible. Increasing value of the a RA appears. Furthermore with increase of a phenomena of RA vanishes and 
again for a approaching to 2 is visible, a = 2 is the well known Gaussian case. The fact that the shape of MFPT 
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Figure 7. Survival probability distribution, G(j,t), (left panel) corresponding to the MFPT(7) curve for a = 2, f3 = 
I, a = l/\/2 with H- = 0,H+ = 8 (cf. Fig. 5) and sample cross-sections (right panel) of G(t) for various 7: '+': 
log 10 7 = —6 (low frequency); 'x': log 10 7 = 0.48 (resonant frequency) and log 10 7 = 10 (high frequency). The results 
were calculated by direct integration of Eq. (1) with the time step dt — 10~ 4 and averaged over N = 10 3 realizations. 
Note for log scale: z-axis (left panel) and y-axis (right panel). The inset in the right panel represent behavior of the 
survival probability, G(t), for small t. 




Figure 8. Survival probability distribution, G(j,t), (left panel) corresponding to the MFPT(7) curve for a — 0.9, /3 = 
1,0'= l/\/2 with H- = 0,H + = 8 (cf. Fig. 4) and sample cross-sections of G(t) (right panel) for various 7: '+': 
l°gio 7 ~ — 6 (low frequency); 'x': log 10 7 w 1 (resonant frequency) and '*': log 10 7 = 10 (high frequency). The results 
were calculated by direct integration of Eq. (1) with the time step dt — 10~ 4 and averaged over N — 10 3 realizations. 
Note for log scale: z-axis (left panel) and y-axis (right panel). The inset in the right panel represent behavior of the 
survival probability, G(t), for small t. 



curve for heavy tailed, skewed distributions (e.g. a = 0.9, (3 = 1) arising from the contact with not equilibrated 
bath is similar to the MFPT curve in the presence of the Gaussian noise (a = 2) arising from the contact with 
equilibrated thermal bath is very interesting and makes the problem of recognition of the underlying noise more 
complicated. It is very intriguing that similar behavior is recovered in two opposite limits. The analysis of MFPT 
only does not allow to recognize the character of underlying noise. The examination of survival probability, G(t), 
or FPTs distribution, /(r), provide a deeper insight into the problem and can be useful in making any statement 
about underlying thermodynamics. However, analysis of G(t) provide more information about the system, it is 
not always sufficient. Therefore further studies are required. 

Like in the Gaussian regime, 3 in the every case under the study the MFPT for small 7 is equal to the average 
value of MFPTs for both barrier configurations, i.e. MFPT(7 -> 0) = i[MFTP(ff_) + MFTP (£/+)]. For large 7 
the MFPT is equal to the MFPT over the average potential barrier, i.e. MFPT(7 -> 00) = MFTP (§(#- + H+)) . 

APPENDIX A. STABLE NOISES AND THEIR REPRESENTATIONS 

Let {Xi,X 2 ,X} are some independent random variables distributed according to the same unknown probability 
distribution. The most intuitive definition 8 of the stable variable, can be written as 

aX 1 + bX 2 = cX + d, (8) 

where = denotes equality in a distribution sense, i.e. the a-stable variables are random variables for which the 
sum of random variables is distributed according to the same distribution as each variable. Real constants c, d 
in Eq. (8) allow for rescaling and shifting of the initial probability distribution. From Eq. (8) the characteristic 
function of a-stable random variables can be determined. 15 The characteristic function of the stable distribution 
can be parameterized in various ways. In the strictly stable parameterization 8, 18 L a ^((;a, /x), a characteristic 
function of the Levy type variables is given by 



4>{k) 



cxp 
exp 



-<j a \k\ a (1 - i/3sign(fc) tan ^) + ifik] , for a ± 1, 



(9) 



-a\k\ (l + i/3-|sign(fc) In |fc|) + ifjbk] , for a = 1, 
with a e (0, 2], (3 € [—1, 1], a e (0, 00), fj, e (—00, 00) and <fr(k) defined in the Fourier space 

<j>{k) = J d(e- lki L a ,p(C;<J,n)- (10) 

Parameter a is called the stability index, (3 describes skewness of the distribution, a is responsible for its scaling 
and fi is a location parameter. Generally, for (3 = fj, = PDFs are symmetric while for (3 = ±1 and a e (0, 1) they 
are totally skewed, i.e. C is always positive only or negative only depending on the sign of skewness parameter /3. 8 
Despite the fact that stable probability distributions asymptotically behave like oc |£|~( a+1 ) analytical expression 
for their probability distributions are known only in the few cases. For a = 2 and any (3 resulting distribution is 
Gaussian 

L 2 , (C; a, it) - — L= cxp (- K-p^T) = AT (/j, 2a 2 ), (11) 



2a^ " V , 

however by the matter of tradition (3 = is usually chosen. For a = 1, (3 = Cauchy distribution is obtained 

^,o(C;^) = ; (c _^ 2 + (t2 , (12) 

and for a = 0.5, (3 = 1 Levy-Smirnoff distribution 

i Vv(C;^M) = (^) i (C-M)-ex P (-^- y ). (13) 

For L a jj((; a, ii) distributions moments of order a exist, i.e. the integral i Qi( g(C; cr, n)( a d( is finite. This 
result in the conclusion that for every a variance of stable distribution does not exist and for a < 1 the average 
value also does not exist. 



Position of the Brownian particle is calculated by direct integration of Eq. (1) with respect to the a-stable 
measure L a jj{s) 



x(t) = - (\v'{x{s))- gr)(s)] ds + f dL a . (j {s) 

J t() J to 

= - f V ± {x(s))ds+ [ dL a . fj {s). 

Jtn Jtn 



(14) 



'to JtQ 

The -L Q ./3 measure in Eq. (14) can be approximated by 8 ' 10 

JV-l 



f(s)dL a ,p(s) 



to 



J2 f(iAs)M a ^ ([iAs, (i + l)As)) 

i=0 
N-l 

/(^A 5 )A S 1/Q Q, 



(15) 



7=0 



where q is distributed according to the L a ^(<;;a, \i = 0), NAs = t — t and M a ^([iAs, (i + l)As)) is the measure 
of the interval [iAs, (i + l)As). 

Stable random variables corresponding to the characteristic function (9), can be generated according to 
the following recipes. 8 ' 18 For a ^ 1 it is necessary to calculate 



with constants B,C,D given by 



sin(a(y + C a ,p)) 



(cos V)- 



cos{V - a{V + C a , p )) 



W 



C a ,(3 



arctan (/3tan(^)) 



a 



D. 



a,f3,a 



= a 



cos ^arctan {ji tan( 



7ra 
~2~ 



For a = 1, <r can be obtained from the formula 



2a 



(-+/3V)tan0O-/?ln 



•^WcosV 

JTW 



(16) 

(17) 
(18) 

(19) 



In the above equations V and W are independent random variables. V is uniformly distributed in the interval 
(—§,§) and W is exponentially distributed with a unit mean. 8,18 

In the simulation, Eq. (14) is numerically examined till x < 1. When x > 1 for the first time the recorded 
value of t = t is an FPT for the problem under the study. From the ensemble of FPTs, r, the MFPT, (r), is 
estimated and the survival probability, G(t), is constructed. Survival probability, G(t), is connected with the 
cumulative distribution of the FPTs 

F(t) = 1 - G(t)| t=T . (20) 
Furthermore, from G(t) the probability density of FPTs, /(t), 7 ' 19 can be calculated 



fir) = -jG(t)\ t = T . 



(21) 
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